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ABSTRACT 
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This  final  report  presents  the  experimental  and  theoretical  work  performed  in  our 
research  program  on  static  and  fatigue  damage  in  high  temperature  composites.  The 
theoretical  part  focused  on  development  and  implementation  of  a  new  Transformation 
Field  Analysis  (TFA)  for  inelastic  laminates,  which  employs  unit-cell,  periodic  array 
models,  as  well  as  averaging  micromechanical  models  together  with  any  constitutive  law  for 
thermo— viscoplastic  and  other  inelastic  deformation.  Moreover,  we  have  initiated 
modeling  of  damage  by  debonding  and  frictional  sliding  at  the  fiber  interface,  using  the 
TFA  computational  scheme  and  the  finite  element  method. 

In  the  experimental  program,  we  examined  the  effect  of  static  and  cyclic  loading  (at 
0.1  and  0.001  Hz)  on  the  inelastic  response  of  and  damage  development  in  (0)4,  (0/90)s  and 
(0/±45/90)s  SiC/Ti  (Sigma/Timetal  21S)  laminates,  and  the  Timetal  21S  matrix,  at  650°C 
and  21°C.  At  the  elevated  temperature,  viscoplastic  deformation  of  the  matrix  can  be 
observed  even  at  relatively  low  applied  stresses.  At  both  temperatures,  reduction  of 
unloading  elastic  moduli  of  the  laminates,  which  indicates  onset  of  damage  by  interface 
decohesion,  also  starts  after  loading  to  relatively  low  stresses.  The  cyclic  loading  rate  has 
no  effect  on  endurance  of  the  (0/9Q)s  and  (0/±45/90)s  laminates  loaded  by  constant 
amplitude  tension  at  650°C.  In  the  (0)4  laminate,  higher  endurance  limits  were  detected  at 
0.1  Hz  than  at  0.001  Hz.  However,  regardless  of  rate  and  layup,  the  total  strain  at  failure 
under  both  static  and  cyclic  loading  at  650°C,  was  measured  at  1±0.1%. 

The  experimental  results  are  interpreted  by  micromechanical  models  developed  in 
the  theoretical  work.  The  effects  of  temperature  and  time  variations  on  the  constitutive 
response  of  the  matrix  and  fiber  are  considered.  The  results  show  that  the  inelastic 
deformation  of  the  titanium  matrix,  assisted  by  extensive  fiber  debonding  in  the  off-axis 
plies,  promotes  stress  transfer  to  the  fibers  in  the  0— degree  plies.  Under  repeated  loads,  the 
fiber  strain  increase  gradually  at  a  certain  rate,  which  depends  on  the  applied  load 
frequency,  and  the  laminate  fails  when  the  fiber  strain  reaches  a  critical  value. 
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1.  INTRODUCTION 
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Fatigue  failure  in  fibrous  composite  laminates  is  typically  associated  with  internal 
stress  redistribution  during  cyclic  loading  that  results  in  stress  buildup  and  overloading  of 
the  fibers  in  certain  plies.  Under  axial  loading  of  laminated  plate  specimens  used  in  most 
experiments,  this  happens  in  the  fibers  of  the  0-degree  plies.  Several  distinct  mechanisms 
may  participate  in  the  stress  redistribution.  For  example,  in  inviscid  elastic— plastic 
systems  at  room  temperature,  matrix  cracking  and/or  fiber— matrix  debonding  in  the 
off-axis  plies  causes  reduction  of  stiffness  of  the  damaged  plies,  so  that  an  increasing  part 
of  the  applied  load  is  supported  by  the  axial  plies. 

A  different  mechanism  operates  in  the  present  SiC/Ti  system  at  elevated 
temperatures,  where  the  matrix  is  viscoplastic  and  thus  its  deformation  depends  on  the 
time  spent  at  a  particular  load  level.  While  both  static  strength  and  isothermal  endurance 
limits  of  such  laminates  are  relatively  high,  it  is  well  known  that  both  the  lack  of 
thermomechanical  compatibility  between  the  fiber  and  matrix,  and  the  relatively  low 
strength  of  the  fiber— matrix  interface,  make  the  laminates  susceptible  to  interface 
decohesion,  fiber  sliding,  matrix  cracking,  and  other  forms  of  distributed  damage. 
Environmental  exposure  is  also  a  factor  at  higher  temperatures.  For  example,  the  hold 
time  at  maximum  stress  in  isothermal  fatigue  tests  performed  in  a  high  temperature  air 
environment  can  severely  decrease  the  fatigue  life  of  a  titanium  matrix  composite  [1].  The 
stress  amplitude,  frequency,  and  the  ratio  of  maximum  to  minimum  stress  have  all  been 
found  to  affect  the  fatigue  behavior  of  unidirectional  composites  under  isothermal  as  well  as 
thermomechanical  loading  conditions  [2-5].  Failure  is  generally  preceded  with  matrix 
transverse  cracking  and  caused  by  overloading  the  fibers. 

Similar  effects  accompanied  with  more  complex  damage  mechanisms  have  been 
observed  in  titanium  matrix  laminates.  Under  isothermal,  and  in— phase  thermomechanical 
fatigue  loading,  debonding  of  the  off-axis  plies  is  the  dominant  damage  mode  [6-8]. 
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Time-dependent  deformation  of  the  matrix  is  significant  under  isothermal  and  in— phase 
loading  conditions  [5,  8—10].  In  contrast,  matrix  flow  is  nearly  absent  under  out— of— phase 
thermomechanical  load  cycles.  In  this  case,  microcracking  of  the  matrix  followed  by  fiber 
breakage  is  the  dominant  damage  mode.  In  any  event,  failure  of  the  composite  is  again 
caused  by  overloading  the  fibers  in  the  0°  plies. 

In  this  report,  we  present  the  results  of  isothermal,  at  21°C  and  650°C,  static  and 
cyclic  loading  experiments,  and  their  theoretical  interpretation,  on  (0)4,  (0/90)s  and 
(0/±45/90)s  Sigma/Timetal  21S  laminates.  Elastic  and  inelastic  response,  damage 
accumulation  detected  by  measuring  changes  in  secant  and  unloading  elastic  moduli,  and 
failure  strains  were  recorded.  While  there  are  many  differences  in  the  behavior  of  the 
tested  laminates,  it  turns  out  that  the  strains  at  failure  are  remarkably  consistent  and 
independent  of  loading  conditions.  At  room  temperature  the  axial  strain  in  the  0°  ply  is  in 
the  range  of  0.5— 0.6%,  while  at  650°C  this  strain,  including  the  thermal  strain,  was  equal 
to  0.9-1. 1%. 

The  experimentally  observed  fatigue  behavior  of  the  Sigma/Timetal  21S  laminates 
have  been  simulated  by  micromechanical  models.  A  summary  of  the  theoretical  approach 
for  analysis  of  fibrous  composite  laminates  is  presented  and  followed  by  numerical 
simulations  of  the  experiments.  The  results  indicate  that  cyclic  loading,  combined  with 
fabrication  and  processing,  produce  complex  internal  stress  states  in  the  plies  and  in  the 
fiber  and  matrix  phases  and  cause  interface  decohesion  by  sliding  or  separation  under 
overall  loads  which  are  only  a  fraction  of  the  ultimate  strength.  Comparison  of  the 
predicted  results  with  experimental  data  suggest  that  fiber  overloading  is  responsible  for 
fatigue  failure  under  constant  temperature.  At  failure,  the  computed  overall  axial  strains 
for  all  laminates  fall  within  a  narrow  band  of  1.0  ±  0.1  %,  which  compares  favorably  with 
the  strain  range  measured  experimentally. 
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2.  EXPERIMENTAL  INVESTIGATION  OF  ISOTHERMAL  STATIC  AND  FATIGUE 
BEHAVIOR  OF  SIGMA/TIMETAL  21S  LAMINATES 

2.1  Materials  and  Specimens 

The  material  studied  herein  is  a  Timetal  21S  matrix  reinforced  by  continuous 
aligned  Sigma  fibers.  The  Timetal  21S  material  is  a  metastable  beta  titanium  alloy  with 
the  composition  (Ti— 15Mo— 3Nb— 3A1— .2Si)  produced  in  foil  form  by  TIMET,  Inc.,  Denver, 
Colorado.  The  Sigma  fiber  is  a  silicon— carbide  monofilament  with  a  structured 
carbon— titanium  boride  coating  manufactured  by  British  Petroleum.  The  fiber  diameter  is 
100/zm  (0.004  in.).  Fabrication  of  silicon-carbide/titanium  composite  panels  typically 
involves  hot  isostatic  pressing  of  alternating  layers  of  titanium  foil  and  tapes  of  aligned 
fibers,  held  in  place  by  molybdenum  wires.  In  particular,  the  Sigma/Timetal  21S 
composite  panels,  supplied  by  Wright  Laboratories,  were  pressed  at  899°C  (1650°F)  under 
a  hydrostatic  pressure  of  103.5  MPa  (15  Ksi)  for  2  hours.  Subsequent  treatment  involved 
aging  at  621°C  (1150°F)  for  8  hours  in  a  vacuum  of  10 '5  torr. 

The  material  was  received  as  flat  plate  panels,  measuring  210  mm  (8.25  in.)  in 
direction  of  the  0°  fiber,  in  three  layups,  (0)4,  (0/90)s,  and  (0/±45/90)s.  C-scan  images  of 
the  laminated  panels  provided  by  the  manufacturer  did  not  reveal  any  major  defects.  The 
fiber  volume  fraction  computed  from  cross  section  images  taken  for  the  received  panels  was 
found  to  be  28.2%  for  the  (0)4  unidirectional  composite,  32.5%  for  the  (0/90)s  cross  ply, 
and  34.1%  for  the  (0/±45/90)s  laminate.  An  unreinforced  plate  manufactured  by  hot 
isostatic  pressing  of  Timetal  21S  foil  was  also  received  for  evaluation  of  the  neat  matrix 
properties.  Specimens  with  straight  edges  were  cut  from  the  panels  using  a  diamond  cutoff 
wheel.  At  3000  rpm  and  a  feed  rate  of  8  mm  per  minute,  this  cutting  method  preserved  the 
laminate  integrity  at  the  cut  section  and  provided  a  smooth  edge.  The  nominal  width  of 
the  specimens  was  10  mm  (0.4  in.)  for  the  neat  matrix  specimens  and  the  (0)4  laminate, 
and  12.5  mm  (0.5  in.)  for  the  (0/90)s  and  the  (0/±45/90)s  laminates.  Surfaces  of  the 
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specimens  were  chemically  cleaned  before  testing.  Aluminum  tabs  were  bonded  to  the  ends 
of  the  specimens  using  3M  single  part  epoxy  adhesive  and  cured  at  107°C  (225°F).  The 
tabs  were  1.6  mm  (1/16  in.)  thick  and  25.4  mm  (1.0  in.)  long,  providing  a  clear  length  of 
160  mm  (6.25  in.)  for  heating  the  gauge  area  and  cooling  the  grips  as  explained  in  the  next 
section. 

2.2  Experiments 

Axial  mechanical  loading  of  the  specimens  was  applied  by  a  MTS  810  uniaxial 
testing  machine  through  screw  loaded  flat  grips.  Heating  of  the  specimens  was  achieved  by 
a  steel  succeptor  tube  which  enclosed  the  specimen  and  was  heated  by  direct  induction 
using  a  5  KW  solid  state  heating  unit.  The  succeptor  tube  was  25.4  mm  (1.0  in.)  in 
diameter  and  100.0  mm  (4.0  in.)  long.  This  allowed  cooling  water  pipes  to  be  installed 
adjacent  to  the  grips.  Temperature  measurements  and  control  were  performed  with 
thermocouples.  Figure  1  shows  the  temperature  profile  measured  along  a  Sigma/Timetal 
specimen  heated  in  the  succeptor  furnace.  A  uniform  temperature  was  maintained  in  the 
50.8  mm  (2.0  in.)  long  center  zone,  while  the  grips  were  kept  at  room  temperature. 

Strain  measurements  were  made  using  a  MTS  high  temperature  extensometer  with 
alumina  rods  25.4  mm  (1.0  in.)  apart.  The  extensometer  rods  have  45°  knife  edges  which 
were  mounted  at  the  specimen  edge  through  holes  drilled  in  the  succeptor  furnace. 
Mounting  the  extensometer  rods  at  the  specimen  edge  rather  than  the  flat  side  prevented 
bending  of  the  specimen  by  the  rods  pressure. 

Both  static  and  fatigue  tests  were  conducted  on  the  neat  matrix  and  the  composite 
laminate  specimens  at  a  constant  temperature  as  described  in  the  subsequent  sections.  The 
load  and  temperature  control,  and  the  data  aquisition  were  performed  by  a  dedicated 
personal  computer. 
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2.3  The  Neat  Matrix 


The  time— dependent  behavior  of  the  Timetal  21S  material  under  static  loading  at 
650°C  was  examined  in  a  multistep  creep  test.  Load  was  applied  at  50  MPa/min  to  stress 
levels  of  50,  75,  100,  125,  and  150  MPa  and  maintained  for  6  hours,  Fig.  2.  The  thermal 
strain  caused  by  heating  the  specimen  from  room  temperature  to  650°C  was  found  to  be 
0.61%.  Under  static  loading  at  650°  the  Young’s  elastic  modulus  was  64.0  GPa  and  the 
yield  strength  30.0  MPa.  The  total  strain  recorded  under  the  multistep  creep  load  is  also 
shown  in  Fig.  2.  The  creep  strain  rates  corresponding  to  each  stress  level  are  shown  in 
Fig.  3.  Data  from  these  and  related  available  experiments  were  utilized  in  finding  matrix 
elastic  and  viscoplastic  properties,  reported  in  Tables  5  and  6. 

Secondary  creep  at  low  stress  levels,  and  tertiary  creep  at  high  stress  levels 
dominate  the  response  to  monotonic  or  sustained  loads  at  650°C.  The  matrix  specimen 
sustained  total  strain  of  10%  in  this  experiment.  Of  course,  in  a  composite  reinforced  by 
stiff  fibers,  the  matrix  strain  will  be  limited  by  the  ultimate  fiber  strain  of  about  1%.  We 
note  that  the  in  situ  matrix  stresses  are  generated  by  a  complex  thermomechancial  loading 
path  induced  by  hot  isostatic  pressing,  cooling  to  room  temperature,  reheating  to  annealing 
temperature,  and  after  another  cooling  cycle,  to  the  test  temperature  of  650°C.  Inelastic 
deformation  during  these  load  cycles,  and  damage  development  that  accompanies  the 
experimental  loading  program,  both  contribute  to  the  final  stress  values.  Figures  2  and  3 
suggest  that  even  at  low  stress  levels,  matrix  time-dependent  behavior  can  be  significant. 

Cyclic  loading  experiments  examined  the  effect  of  loading  frequency  and  stress  range 
on  strain  response  of  the  Timetal  21S.  Figure  4  shows  the  stress  history  applied  to  a  neat 
Timetal  21S  specimen  at  650°C.  Axial  stress  cycles  with  amplitude  Smax  =  40,  100,  200, 
and  250  were  applied,  while  the  stress  ratio  R  =  Smax/Smin  was  kept  constant  at  0.1.  The 
loading  frequency  was  changed  from  fast  at  1.0  Hz  to  slow  at  0.01  Hz  as  indicated  in  Fig.  4. 
Figure  6  shows  the  total  strain  measured  during  the  entire  history  applied  in  this 
experiment.  The  Timetal  matrix  exhibits  significant  ratcheting  under  cyclic  loading.  As 
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in  the  static  loading  case,  the  cyclic  creep  strain  response  shown  in  Fig.  5  tends  to  stabilize 
if  the  maximum  stress  does  not  exceed  200  MPa.  Tertiary  cyclic  creep  is  observed  under 
load  cycles  with  the  stress  amplitude  of  250  MPa.  As  indicated  in  Fig.  6,  where  the  strain 
developed  under  each  stress  range  is  shown,  the  ratcheting  rate  of  the  Timetal  matrix  is 
affected  by  the  stress  amplitude  and  loading  frequency. 

2.4  Laminates 

Three  laminate  layups,  (0)4,  (0/90)s,  and  (0/±45/90)s,  were  tested  under  static  and 
cyclic  loads  applied  at  constant  temperatures  of  650°C  and  21°C.  The  elevated 
temperature  experiments  were  designed  to  examine  the  effect  of  the  time— dependent 
behavior  of  the  matrix  and  of  distributed  damage  on  the  overall  response.  The 
room— temperature  tests  were  used  only  to  evaluate  damage-induced  changes  of  elastic 
moduli  in  virgin  laminate  specimens  loaded  to  different  stress  levels,  and  of  specimens  that 
survived  certain  cyclic  loading  programs  at  650°C. 

Static  tensile  tests  at  650°C  are  summarized  in  Table  1,  and  plotted  in  Fig.  7.  Two 
loading  rates,  0.278  MPa/s  and  27.8  MPa/s  were  employed.  Stress  rates  of  the  same  order 
were  applied  in  the  cyclic  loading  tests  described  subsequently.  The  elastic  modulus, 
ultimate  strength,  and  maximum  strain  at  failure  are  indicated  in  Table  1  together  with 
the  thermal  strains  caused  by  heating  to  650°C.  The  measured  stress— strain  response  for 
the  three  laminates  is  shown  in  Fig.  7.  Constitutive  behavior  of  the  Sigma/Timetal  2 IS 
composite  at  650°C  is  clearly  rate— dependent.  Some  rate  effects  in  the  elastic  modulus  and 
in  the  maximum  stress  and  strain  are  observed,  particularly  in  the  (0/90)s  and  (0/±45/90)s 
laminates,  Table  1.  In  general,  slow  loading  leads  to  lower  matrix  stresses  and  higher  fiber 
stresses,  and  thus  to  lower  overall  strength,  as  seen  in  Table  1.  However,  regardless  of  the 
loading  rate,  the  combined  mechanical  and  thermal  maximum  strain  sustained  by  the 
laminates  was  equal  to  1.0  —  1.2%. 
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The  laminate  thermal  strains  at  the  test  temperature  are  also  shown  in  Table  1; 
they  are  25  to  30%  lower  that  the  matrix  thermal  strain  of  0.61%.  Of  course,  this 
reduction  is  caused  by  the  constraint  imposed  by  the  fibers,  particularly  by  the  0— degree 
fibers  which  thus  support  a  residual  stress.  For  example,  the  computed  (Section  3.3)  fiber 
axial  thermal  stress  after  heating  of  initially  stress— free  Sigma/Timetal  21S  laminates  to 
650°C  is  equal  to  634  MPa  in  the  unidirectional  laminate;  in  the  (0/90)s  and  (0/±45/90)s 
la.minat.ps,  the  computed  0— degree  ply  fiber  axial  stress  is  equal  to  800  MPa.  However, 
these  stresses  are  superimposed  with  the  thermal  residual  stresses  generated  during 
fabrication.  Cooldown  from  processing  temperature  causes  high  compressive  axial  residual 
stress  in  the  fiber,  hence  the  net  axial  stress  in  the  fiber  after  reheating  to  650°C  is 
compressive.  For  the  Sigma/Timetal  21S  system,  this  stress  was  computed  as  equal  to 
—136  MPa  in  the  unidirectional  composite  (Table  7),  and  —55  MPa  in  either  the  (0/90)s  or 
the  (0/±45/90)s  laminate  (Table  8). 

The  effect  of  the  time-dependent  response  of  the  matrix  on  the  overall  response  of 
Sigma/Timetal  21S  laminates  was  examined  by  multistep  creep  tests.  Each  specimen  of 
the  (0)4,  (0/90)s  and  (0/±45/90)s  laminates  was  heated  to  650°C,  and  then  loaded  axially 
at  50  MPa/min,  to  a  specific  tensile  stress  magnitude  which  was  sustained  for  24  hours. 
Figure  8  shows  the  overall  axial  strain  recorded  in  a  multistep  creep  test  for  two  (0)4 
unidirectional  specimens,  one  with  sustained  axial  stress  of  50,  100,  150,  and  200  MPa,  the 
other  with  a  stress  of  300,  400,  and  500  MPa.  A  limited  amount  of  overall  creep  strain  was 
observed  in  this  composite  system  due  to  the  axial  constraint  imposed  on  the  matrix  by  the 
elastic  Sigma  fibers.  In  elastic  deformation  at  650°C,  the  ratio  of  the  matrix  to  overall 
axial  stress  is  equal  to  about  0.4.  Under  viscoplastic  deformation,  the  matrix  stresses  are 
gradually  reduced  by  creep  strains  at  decreasing  rates.  Figure  8  shows  that  the  second 
specimen  failed  at  1%  overall  strain,  after  a  short  creep  period  at  500  MPa.  The  total 
creep  strain  accumulated  in  this  specimen  was  computed  as  equal  to  0.2%. 
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Higher  creep  strains  were  found  in  Fig.  9  for  the  (0/90)s,  and  in  Fig.  10  for  the 
(0/±45/90)s  laminate,  due  to  the  presence  of  the  more  compliant  off-axis  plies.  However, 
as  in  the  unidirectional  composite,  failure  in  the  (0/±45/90)s  laminate  occurred  at 
approximately  1%  overall  axial  strain.  We  note  that  the  creep  test  of  the  (0/±45/90)s, 
Fig.  10,  was  interrupted  after  about  20  hours  (point  "A"  in  Fig.  10)  due  to  power  failure, 
and  was  resumed  at  the  sustained  stress  of  100  MPa.  This  caused  the  discontinuity  in  the 
strain  rate  found  at  point  "A"  in  Fig.  10. 

Cyclic  tensile  tests  were  conducted  on  the  three  laminates  at  various  maximum 
stresses.  The  ratio  of  the  maximum  to  minimum  applied  stresses,  R  =  Smax/Smin,  was  kept 
constant  at  0.1.  Two  load  frequencies  were  selected,  a  fast  cycle  at  0.1  Hz,  and  a  slow  cycle 
at  0.001  Hz.  Samples  of  the  raw  experimental  data  obtained  for  laminates  are  plotted  in 
Figs.  11—22.  In  a  typical  fatigue  test,  a  specimen  was  first  heated  to  650°C  in  6  minutes 
and  kept  at  this  temperature  for  4  minutes.  The  specimen  was  then  loaded  axially  to  the 
mean  stress  in  one  minute.  Next,  the  mean  stress  was  kept  constant  for  one  minute  before 
application  of  the  load  cycles. 

The  temperature  and  stress  histories  for  the  unidirectional  (0)4  laminate  under  the 
slow  0.001  Hz,  and  the  fast  0.1  Hz  cycles  are  shown  in  the  insets  of  Figs.  11  and  12  which 
display  the  corresponding  strain  histories.  The  stress  amplitudes  were  450  MPa  and 
850  MPa,  respectively.  The  recorded  stress— strain  response  is  shown  in  Figs.  13  and  14. 
Ratcheting  of  the  titanium  matrix  is  reflected  in  the  overall  response.  Under  slow  load 
cycles,  higher  viscoplastic  strains  are  generated  by  each  load  cycle,  Fig.  11,  hence  in  only 
few  cycles,  the  computed  in  situ  matrix  stress  state  becomes  nearly  quasi— isotropic,  and 
the  overall  response  stabilizes.  In  contrast,  several  hundred  cycles  are  required  to  reach  a 
stable  overall  response  when  fast  load  cycles  are  applied.  Similar  behavior  was  found  in  the 
(0/90)s  cross  ply,  Figs.  15—18,  and  in  the  (0/±45/90)s  laminate,  Figs.  19—22. 

While  in  all  the  tested  samples  the  overall  strain  at  failure  is  equal  to  about  1%,  the 
strain  accumulation  per  cycle  depends  on  the  loading  frequency.  This  is  particularly 
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evident  in  the  (0/±45/90)s  sample  results  shown  in  Figs.  19  and  20  where  the  amplitude  of 
250  MPa  is  constant  while  the  loading  frequency  is  the  only  variable.  More  than  2800 
cycles  at  0.1  Hz  were  required  to  bring  the  overall  strain  to  the  ultimate  value,  Fig.  20. 
Under  the  same  overall  stress  applied  at  0.001  Hz,  the  specimen  failed  in  only  215  cycles. 

Fatigue  experiments  conducted  on  the  Sigma/Timetal  21S  laminates  are  summarized 
in  Table  2.  The  S— N  data  appear  in  Fig.  23,  which  also  shows  data  points  for  failure  in  the 
first  cycle,  that  were  not  included  in  Table  2.  The  observed  endurance  limits  for  (0/90)s 
and  (0/±45/90)s  laminates  are  similar.  The  (0/90)s  crossply  shows  a  somewhat  longer 
fatigue  life  than  the  (0/±45/90)s  laminate.  Although  the  S-N  data  under  the  fast  and  slow 
load  cycles  ai  similar,  the  fatigue  life  of  either  laminate  is  longer  at  the  higher  loading 
frequency.  As  discussed  earlier,  this  is  caused  by  the  slower  growth  of  the  matrix  creep 
strains  under  fast  load  cycles.  According  to  the  data  in  Fig.  23,  the  endurance  limit  of  the 
(0/90)s  crossply  is  300  MPa;  for  the  (0/±45/90)s  laminate  it  is  200  MPa. 

The  unidirectionally  reinforced  specimens  show  a  more  pronounced  effect  of  the 
loading  frequency  on  the  fatigue  life,  this  was  also  reported  by  Nicholas  and  Russ  [1]. 
However,  the  S-N  data  in  Fig.  23  also  indicate  significant  differences  between  fatigue  lives 
of  two  batches  of  nominally  similar  unidirectional  laminates.  Specimens  cut  from  the 
B— piite  had  an  order  of  magnitude  longer  fatigue  life  than  the  specimens  cut  from  the 
A— plate. 

Damage  development  was  monitored  by  measuring  the  unloading  elastic  moduli  and 
secant  moduli  of  selected  specimen.  Table  3  shows  room— temperature  results  obtained  for 
the  (0/90)s  layup.  One  virgin  or  reference  specimen  was  loaded  at  a  slow  rate  to  the  stress 
levels  indicated  in  the  first  column  of  Table  3.  At  each  level,  the  initial  unloading  elastic 
modulus  was  measured  during  a  complete  unloading  cycle.  This  was  followed  by  loading  to 
a  higher  level  and  another  unloading.  The  Young’s  modulus  measured  during  first  loading 
was  equal  to  171  GPa.  Table  3  also  shows  room-temperature  unloading  moduli  measured 
by  the  said  procedure  on  two  specimens  that  survived  110000  and  113750  cycles  at  650°, 
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respectively.  Ultimate  failure  stress  levels  for  the  three  specimens  are  also  given  in  the 
table.  Table  4  presents  data  from  similar  tests  performed  on  (0/±45/90)s  laminates,  where 
the  Young’s  modulus  measured  during  first  loading  was  equal  to  178  GPa. 

In  both  laminates,  the  unloading  moduli  are  lower  than  the  respective  initial 
(loading)  magnitudes,  and  decrease  with  the  rising  applied  stress.  This  is  indicative  of 
damage,  such  as  interface  debonding  and  sliding  in  the  virgin  reference  samples,  and  also  of 
matrix  and  fiber  cracking  in  those  exposed  to  isothermal  fatigue,  and  of  progressive  opening 
of  new  or  pre-existing  separations  with  higher  applied  stresses.  Indeed,  matrix  cracks  were 
detected  on  micrographs  of  sectioned  samples,  and  so  was  interface  decohesion,  but  no 
attempt  was  made  to  quantify  such  observations.  It  is  clear,  however,  that  the  present 
laminates  are  very  prone  to  damage,  even  at  static  loads  well  below  ultimate  strength. 

Damage  development  was  also  detected  by  measuring  changes  in  the  secant  moduli 
computed  from  the  hysteresis  loops,  such  as  those  appearing  in  Fig.  13.  Results  obtained 
on  two  (0)4  specimens  from  plate  "A"  cycled  at  0.001  Hz  are  shown  in  Fig.  24,  where  the 
secant  modulus  reduction  is  related  to  the  overall  axial  strain.  The  initial  10-15% 
reduction  of  the  secant  modulus  was  probably  caused  by  viscoplastic  deformation  of  the 
matrix,  which  stabilizes  after  a  few  cycles  at  about  0.85—0.9%  overall  axial  strain. 
Additional  cycles  caused  further  reductions  in  the  secant  modulus  that  continue  until 
failure;  this  indicated  evolution  of  distributed  damage.  In  contrast,  the  reduced  secant 
modulus  of  a  (0)4  specimen  from  plate  "B",  which  survived  1750  cycles  at  a  higher 
maximum  stress  of  600  MPa,  remained  stable  up  to  0.96%  overall  axial  strain. 

Total  axial  strain  prior  to  failure  was  measured  on  all  fatigue  specimens;  for  all  three 
laminates  the  magnitudes  remained  within  the  range  of  1.0  ±  0.1  %,  Fig.  25.  Since  the 
residual  thermal  and  inelastic  strains  are  mostly  relieved  by  damage,  this  suggests  that 
failure  of  the  tested  composites  was  caused  by  overloading  the  Sigma  fibers  in  the  0°  ply 
which  can  sustain  an  ultimate  strain  of  approximately  1%  at  650°C.  This  critical  strain 
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appears  to  be  a  useful  design  criterion,  providing  that  the  contribution  due  to  damage  can 
be  properly  accounted  for. 
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3.  MODELING  OF  THERMOMECHANICAL  BEHAVIOR  OF  HIGH  TEMPERATURE 

METAL  MATRIX  COMPOSITE  LAMINATES 


3.1  Constitutive  Equations  of  the  Fiber  and  Matrix 

Both  constituents  are  regarded  as  homogeneous  solids  with  known, 
temperature-dependent  thermo-elastic  properties;  both  phases  are  isotropic  in  the  present 
application,  but  the  analysis  can  accommodate  any  actual  material  symmetry.  The  fiber 
remains  elastic  until  failure.  The  matrix  is  elastic  if  loaded  within  a  certain  equilibrium 
surface  in  the  stress  space,  but  becomes  viscoplastic  if  loaded  beyond  this  surface.  The 
inelastic  response  of  the  matrix  is  derived  here  from  a  unified  viscoplastic  theory  described 
by  Bahei— El— Din  et  al.  [11]  and  Shah  [12].  In  either  phase,  the  total  strain  rate  is 
additively  decomposed  into  elastic,  thermal,  and  inelastic  components, 


cij  = 


•  e  .  t  .in 

+  €::  +  €z 


'1J 


“1J 


“1J 


The  elastic  and  thermal  rates  are  found  from 
*ij  =  Mijki(0)  <^ki  > 

til  =  KdMyuM/dtf)  <rkl  +  fy  M)  9 


(1) 


(2) 

(3) 


0 

where  9  is  current  temperature,  Mijkl(0)  is  the  elastic  compliance,  is  the  Kronecker 
tensor  and  (3(6)  is  linear  thermal  expansion  coefficient;  (3)  follows  from  (A— 1)  in  Appendix. 

The  inelastic  strain  is  found  as  a  function  of  the  overstress  measured  from  an 
equilibrium  yield  surface  which  contains  all  stress  states  that  can  be  reached  from  the 
current  state  by  purely  elastic  deformation  [11,  12].  The  equilibrium  surface  is  given  by 

«  =  |  (slj  -  »y)  <sh  -  au)  -  <Y  +  Q)J  =  °,  (4) 
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4c 

where  Sy  is  the  deviatoric  equilibrium  stress  tensor,  ay  denotes  the  center  of  the  yield 
surface,  Y  =  Y(0)  is  the  initial  yield  stress  in  tension,  independent  of  the  loading  rate,  and 
Q  =  Q(0)  is  an  isotropic  hardening  function. 

For  a  given  stress  state  Sy  outside  the  surface  (4),  there  exists  an  equilibrium  stress 

4c  4c 

Sjj  which  satisfies  (4)  such  that  the  stress  points  Sy  and  Sy,  and  the  center  of  the 
equilibrium  surface  ay  are  collinear; 


4c 


2(Y+  Q)2 

3(skr  ak*)(skr  ak*) 


i 

(Sy  -  ay)  +  ay  . 


(5) 


The  effective  overstress  R  is  a  measure  of  the  distance  between  the  actual  stress  point  Sy 
and  the  equilibrium  stress  point  Sy.  It  vanishes  if  the  stress  point  lies  on,  or  falls  inside 
the  yield  surface.  In  particular, 

R  =  [|  (Sy  -  Sy)(Sy  -  Sy)]*  if  g(Sy  -  ay)  >  0  , 

(6) 

R  =  0  if  g(Sy  -  ay)  <  0  . 


The  inelastic  strain  rate  is  then  written  as 


£§  =  ^3/2)k(«)RP'9’nij> 


(?) 


where  the  functions  k(0)  and  p(0)  are  material  parameters  and  ny  is  a  unit  normal  to  the 
yield  surface  (4)  at  the  current  equilibrium  stress  point.  From  (4)  one  finds  this  as 


nij 


(■g  ~  aij) 

[(sk*“  ak()  (sk{~  ak/)]^ 


(sij  ~  alj) 
(Y  +  Q) 


(8) 


The  evolution  equation  for  Q  which  includes  the  effect  of  inelastic  deformation  and 
thermal  recovery  on  the  yield  stress  is  given  as 
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Q  =  q(*)[<Wf>)-C!]71"-br(<*)  |Q-Q,(<0|<V”'U[Q-Qr(«)].  (9) 

The  functions  q( 9),  Qa(0),  br( 0),  Qr( 0),  and  nr(0)  are  material  parameters,  total  (Qr  =  0) 

or  partial  (Qr  ^  0)  thermal  recovery  is  represented  by  the  second  term  in  (9),  and  7ln  is  the 
effective  inelastic  strain  rate; 
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(10) 


In  analogy  to  (9),  and  permitting  complete  thermal  recovery  of  kinematic 
hardening,  the  evolution  equation  for  the  center  ay  of  the  yield  surface  can  be  written  as 

aij  =  M  vij  ~  M#)  a  ay  ,  a  =  (ak/  ak/)*  ,  (11) 


where  vr(  0)  and  wr(0)  are  material  parameters.  The  unit  tensor  */y  defines  the  direction  of 
translation  of  the  yield  surface  in  the  deviatoric  stress  space,  and  can  be  specified  according 
to  the  hardening  rules  applied  in  rate-independent  plasticity  theories.  If  the  Phillips 
hardening  rule  is  selected,  then 


"ij  =  Sij/(sk/sk^ 


VU  =  nii 


if  Sy  ^  0  , 


if  Sy  =  0 


(12) 

(13) 


The  factor  /i  in  (11)  is  found  from  Prager’s  consistency  condition  g  =  0,  with 
translation  of  the  yield  surface  specified  by  the  first  term  in  (11).  The  result  is 


k  =  A W  k(«)  RP<9’  [H(0)  -  q(0)[Q,(0)  -  Q]\/ukl  ^ .  (14) 

A  variant  of  the  two-surface  plasticity  theory  by  Dafalias  and  Popov  [13]  is  used  to 
describe  evolution  of  the  instantaneous  tangent  modulus  H: 

H(0)  =  Ho(0)  +  h(0)  [{/(«,  -S)f  ,  (15) 

{  =  [I  (Sij  ~  sij)(s«  -  8tj)]  ^  >  <16) 


17 


* 

where  6  is  the  distance  between  the  stress  point  s^  on  the  equilibrium  surface  and  a  point 
sy  on  the  bounding  surface  with  normal  equals  n^,  and  S0  is  the  value  of  6  at  the  onset  of 
inelastic  deformation.  When  the  equilibrium  stress  point  lies  on  the  bounding  surface,  the 
plastic  tangent  modulus  H(0)  assumes  the  asymptotic  value  Ho(0),  which  together  with  the 
parameters  h(0)  and  m(0)  need  to  be  determined  experimentally.  In  analogy  with  the 
equilibrium  yield  surface,  thermal  recovery  of  isotropic  as  well  as  kinematic  hardening  of 

t 

the  bounding  surface  can  be  included  in  the  model.  This  is  omitted  here  for  brevity.  We 
only  mention  that  the  recovery  terms  for  isotropic  and  kinematic  hardening  of  the 
bounding  surface  assume  a  form  similar  to  those  suggested  above  for  the  yield  surface,  but 
with  new  material  parameters  which  we  denote  by  script  letter.  For  example,  the  Mises 
form  of  the  bounding  surface  is  written  as 

/  =  |(sij  -  +  •2)2  =  0  »  (17) 

where  is  the  center  of  the  bounding  surface  in  the  deviatoric  stress  space,  /=  /(0)  is 
the  bounding  stress  given  by  the  intersection  of  the  asymptotic  part  of  the  uniaxial 
stress— plastic  strain  curve  and  the  stress  axis,  and  2  =  £  (0)  is  an  isotropic  hardening 
function. 

In  the  Sigma/Timetal  21S  system,  the  fiber  was  regarded  as  elastic  and  the  matrix 
as  thermo— viscoplastic.  Elastic  moduli  of  the  Sigma  fiber  were  assumed  to  be  the  same  as 
those  of  the  SCS6  fiber,  and  the  material  parameters  for  the  Timetal  21S  were  derived,  in 
part,  from  experiments  performed  at  NASA  Langley  Research  Center  [14,  15].  Isotropic 
hardening  of  the  equilibrium  (4)  and  bounding  (17)  surfaces  was  neglected;  thermal 
recovery  in  (9)  and  (11)  was  minor.  The  material  parameters  remaining  after  reductions 
appear  in  Tables  5  and  6. 
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3.2  Micromechanical  Models  of  Composite  Laminates 

3.2.1  Model  selection 

Among  the  available  approaches  to  modeling  of  inelastic  composite  materials,  we 
use  here  the  transformation  field  analysis  (TFA)  of  the  inelastic  lamina  [16]  in  conjunction 
with  the  laminated  plate  theory  for  symmetric  laminates  subjected  to  uniform  in— plane 
overall  stresses,  out-of-plane  normal  stress,  and  changes  in  temperature.  When  used  with 
a  selected  micromechanical  model,  the  TFA  method  provides  piecewise  uniform 
approximations  of  the  phase  stress  and  strain  fields  that  are  then  utilized  in  finding  the 
overall  instantaneous  stiffness,  compliance,  and  the  thermal  stress  and  strain  tensors  of  the 
plies  and  laminate.  The  governing  equations  of  the  TFA  method  are  described  in  some 
detail  below.  These  equations  were  incorporated  in  a  finite  element  procedure  [17],  as 
summarized  in  the  sequel,  and  implemented  in  the  VISCOPAC  program  [18]  for  a  class  of 
averaging  micromechanical  models.  Specifically,  the  vanishing  fiber  diameter  [19,  20],  the 
self-consistent  [21],  and  the  Mori-Tanaka  [22,  23]  models  are  included.  Work  in  progress 
uses  the  Periodic  Hexagonal  Array  (PHA)  model  [24,  25]  that  provides  piecewise  uniform 
approximations  of  the  local  fields  in  the  matrix,  fiber,  and  at  their  interfaces,  and  can  thus 
follow  more  closely  their  evolution  along  the  local  loading  path. 

3.2.2  Transformation  field  analysis 

A  representative  volume  of  a  fibrous  composite  ply  is  considered,  under  uniform 
overall  stress  or  strain,  and  a  uniform  change  in  temperature.  The  goal  is  to  estimate  the 
stress  and  (small)  strain  averages  caused  in  the  elastic  fiber  and  inelastic  matrix  by  the 
applied  loads.  In  the  transformation  field  analysis,  these  averages  are  found  by 
superposition  of  two  load  systems,  one  consisting  of  the  applied  mechanical  loads,  in  the 
form  of  uniform  overall  stress  or  strain,  and  one  associated  with  internal  eigenstrains  that 
include  thermal  strains,  any  inelastic  strains,  and  also  equivalent  eigenstrains  designed  to 
simulate  damage  development,  e.  g.,  by  fiber  sliding  and  debonding.  This  identification  of 
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thermal,  inelastic  and  damage— induced  deformations  with  eigenstrains,  implies  that  both 
load  systems  act  on  an  elastic,  undamaged  composite.  Of  course,  except  for  the  thermal 
strains,  the  eigenstrain  load  system  is  not  known  a  priori,  it  is  evaluated  incrementally 
from  certain  governing  equations  provided  by  the  transformation  analysis  method. 

The  procedure  has  been  described  in  several  references,  for  example  in  [16]  and  [26]. 
An  exposition  of  these  results,  that  is  particularly  easy  to  implement,  uses  symbolic 
notation  where  boldface  lower  case  letters  represent  symmetric  second-order  tensors 
written  as  (6x1)  matrices,  and  boldface  capital  letters  denote  symmetric  fourth-order 
tensors  written  as  (6x6)  matrices.  Connections  with  the  notation  used  in  Section  2  are 
easily  established.  For  example,  the  strain  tensor  e^,  with  the  symmetry  =  e^,  is 
written  as  e  =  {en  e22  e33  2e23  2e31  2e12}T;  fourth-order  tensors  having  at  least  the 
symmetries  AAjkl  =  Aj ikl  =  Aijlk  are  reduced  to  (6x6)  matrices  A,  such  that 
AA  *=  A  *A  =  I,  the  identity  matrix. 

The  selected  representative  volume  V  of  a  composite  ply  is  assumed  to  contain 
many  fibers  and  be  loaded  by  an  overall  stress  a  or  strain  e,  both  macroscopically  uniform. 
A  certain  subdivision  of  V  into  subvolumes  Vr,  Vs  is  chosen  (r,  s  =  1,  2,  ...  N).  In  what 
follows,  N  =  2,  and  Vr,  Vs  will  be  either  fiber  or  matrix  volumes  Vf  or  Vm,  (r,  s  =  f,  m). 

The  elastic  constitutive  relations  for  the  fiber  and  matrix  phases  with  initial  strains 
or  eigenstrains  /i r  and  eigenstresses  Ar  =  —  Lr  are  written  as 

aT  =  Lr  cr  +  Ar  eT  =  Mr  <rr  +  ^ .  (18) 

The  Lr  and  Mr  =  l/  are  phase  elastic  stiffness  and  compliance  matrices,  discussed  in  the 
Appendix.  The  is  the  total  eigenstrain  that  would  remain  in  the  phase  volume  Vr  after 
complete  local  unloading  to  zero  stress,  whereas  the  eigenstress  Ar  is  seen  to  represent  the 
total  stress  caused  in  a  fully  constrained  volume  Vr  by  the  eigenstrain  fi r.  It  then  follows 
that  Ar  =  -  Lr  .  Since  the  Lr  and  Mr  depend  only  on  elastic  moduli,  assumed  to  be 
constant  at  a  given  temperature,  (18)  can  be  readily  converted  into  a  rate  form, 
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®r  =  Lr  «r  +  Ar 


cr  =  M  r  ar  +  . 


(19) 


With  reference  to  (1),  (3)  and  (7),  the  eigenstrain  rate  that  represents  the  thermal  and 
inelastic  deformation,  may  be  written  as 

•  •  c  i  •  •  t  •  in  .  > 

e=e+/i  fi=  e  +  €  .  (20) 

The  superposition  of  the  two  load  states  employed  in  the  transformation  field 
method  provides  the  following  form  of  the  local  fields, 

=  e  ^  ®rs  As  > 

s 

r,  s  =  f,  m  (21) 

K  =  Br  o  +  X  ' 

S 

The  e  and  a  denote  the  uniform  overall  strain  and  stress  rates  applied  to  the  representative 
volume,  and  A,,  and  Br  are  mechanical  strain  and  stress  concentration  factor  matrices.  The 
Drs  and  Frs  are  transformation  concentration  factors,  which  evaluate  the  contribution  to 
the  local  field  in  Vr,  caused  by  a  uniform  transformation  field  of  unit  magnitude  present  in 
Vs.  Both  the  mechanical  and  transformation  concentration  factors  depend  only  on  elastic 
moduli  and  phase  geometry,  and  thus  remain  constant,  except  under  large  changes  in 
geometry  or  temperature  which  affect  the  elastic  moduli.  In  the  latter  case,  the 
concentration  factors  are  recalculated  to  reflect  such  changes,  c.  f.,  §3.2.3  below. 

Of  course,  (19)  can  be  applied  to  any  number  of  subdivisions  of  the  representative 
volume.  However,  we  limit  our  attention  to  phase  average  fields  in  two— phase  solids. 
Then,  the  concentration  factors  can  be  found  in  a  very  simple  form  [16,  27,  28],  involving 
only  overall  elastic  moduli  L,  or  compliances  M,  and  the  known  phase  elastic  moduli  and 
volume  fractions, 
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A,  =  (Lr  —  LS)"‘(L  —  Ls)/cr , 

r,  s  =  f,  m  (22) 

Br  =  (Mr  —  MS)"‘(M  —  Mg)/cr  , 

D,f  =  (I  -  4)  (L,  -  LJ'1  Lt  Drm  =  -  (I  -  Ar)  (Lf  -  L.)"  Lm  , 

(23) 

Frf  =  (I  -  Br)  (M,  -  Mm)->  Mf  F„  =  -(I- B,)  (Mf  -  MJ*'  Mu  . 

The  overall  elastic  stiffness  can  be  either  measured,  bracketed  by  the  Hashin— Rosen  or 
other  exact  bounds,  or  estimated  by  several  approximate  methods;  in  this  work  we  use  the 
Mori— Tanaka  estimates  discussed  in  the  Appendix. 

The  governing  equations  for  the  unknown  local  eigenstrain  rates  can  be  now 
completed  by  first  substituting  for  the  concentration  factors  from  (22)  and  (23),  and  for  the 
eigenstrains  from  (20).  In  particular,  the  /t  in  (20)  is  obtained  from  (3)  and  (7),  converted 

into  an  eigenstress  using  the  relation  A  =  —  L  ji  that  follows  from  (18)  or  (19),  and 
substituted  into  (21 2).  Therefore,  the  eigenstrain  rates  are  replaced  by  expressions 
involving  increments  of  temperature  and  of  local  stress  rates,  and  (212)  is  converted  into  a 
system  of  differential  equations  for  evaluation  of  the  local  stresses.  In  a  similar  way,  (21 1) 
can  be  converted  into  a  system  of  equations  for  the  local  strains.  Solution  techniques  for 
either  system  are  readily  available  [26]. 

Once  the  averages  of  the  local  stress  or  strain  rates  are  known  at  a  particular  point 
of  the  loading  path,  the  inelastic  strains  and  other  local  eigenstrains  can  be  found  from  the 
phase  constitutive  relations.  Overall  averages  of  the  total  local  stress  and  strain  field  rates 
follow  from  the  connections 

°  =  X Cf  ^ 


6  = 


c„  e 


r  > 


(24) 
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(25) 


while  the  overall  transformation  fields  are  found  from, 


A  =  £  cr  A?  Ar 
r 


cr  Br  ^ 


The  overall  constitutive  relations  of  the  ply  are  then  written  in  a  form  similar  to 
(18)  and  (19), 


o 

a  =  L  e  +  A 


T  *0  > 

a  =  L  e  +  A 


o 

£  =  M  a  +  fi 

(26) 

■  _  _  •  0  • 
e  =  M  a  +  fi , 


where  cr  or  e  are  the  overall  applied  rates,  the  overall  elastic  stiffness  L  and  compliance  M 
are  obtained  as  in  (22),  from  experimental  measurements,  bounds,  or  estimates. 

The  transformation  analysis  procedure  can  be  used  with  any  physically  acceptable 
phase  constitutive  relations  and  micromechanical  technique  for  elastic  composites, 
including  the  unit-cell  models  [24,  25],  and  also  in  laminate  analysis.  The  only  difference 
resulting  from  the  various  model  selection  is  in  the  number  of  subdivisions  of  the 
representative  volume  (r  =  1,  2,...N),  and  in  evaluation  of  the  concentration  factor  tensors 
for  each  subvolume  Vr  .  Of  course,  if  more  local  volumes  are  used,  the  piecewise  uniform 
approximation  of  the  local  fields  found  from  (21)  should  be  more  responsive  to  the  inelastic 
constitutive  relations  of  the  phases  and  thus  more  accurate.  However,  for  given 
constitutive  relations,  the  structure  of  the  governing  equations  and  the  solution  procedures 
are  identical  for  all  models. 

This  approach  differs  from  that  where  the  inclusion  problem  is  solved  for  an  elastic 
matrix  whose  compliance  is  identified  with  the  instantaneous  viscoplastic  compliance 
derived  from  the  average  matrix  stress  history  [29].  The  latter  calls  for  repetitive  solutions 
of  the  inelastic  inclusion  problem  in  a  generally  anisotropic  matrix,  not  needed  in  our 
approach;  it  also  may  violate  (25).  On  the  other  hand,  if  the  concentration  factors  in  (21) 
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are  found  using  a  finite  element  analysis  of  a  subdivided  unit  cell  model,  such  as  the  PHA 
model  [24,  25],  then  the  solution  of  (21)  coincides  with  an  inelastic  finite  element 
solution  [26]. 

3.2.3  Modeling  of  interface  damage 

Fiber  debonding  and/or  sliding  are  usually  the  initial  damage  states  in  fibrous 

composite  plies  subjected  to  off-axis  loading  of  sufficient  magnitude.  The  interface 

strengths  are  typically  quite  small.  For  example,  Clyne  and  Withers  [30]  give  interface 

shear  strength  magnitudes  in  the  range  of  50-100  MPa;  normal  strength  values  do  not 

appear  to  be  available,  but  are  expected  to  be  lower.  Cyclic  shear  sliding  tends  to  reduce 

the  number  of  adhering  asperities  and  thus  the  apparent  interface  shear  strength. 

Therefore,  transverse  tension  and  shear,  and  longitudinal  shear  stresses  in  plies  can  initiate 

interface  damage.  Of  course,  processing  and  heat  treatment  may  affect  the  actual 

magnitudes;  for  example,  growth  of  brittle  interfacial  reaction  layers  may  promote  early 

debonding  and  matrix  crack  growth,  whereas  treatments  that  promote  favorable  residual 

stress  states  may  increase  the  apparent  strength. 

However,  as  long  as  some  adhering  asperities  are  left  at  the  interface  after 

debonding,  they  should  provide  for  shear  stress  transfer  that  may  prevent  fiber  pullout. 

For  example,  the  average  surface  shear  (r)  and  normal  ( a )  stresses  in  the  fiber  are 

2 

connected  by  the  force  equilibrium  relation  2toLt  =  7rr  a  over  a  length  L  of  a  fiber  with 
radius  r,  i.  e.,  <r/r  =  2L/r.  If  L  =  25  mm  is  selected  as  a  transfer  length,  and  the  radius  of 
the  Sigma  fiber  r  =  50  /on,  then  one  finds  that  a  =  1000  r.  This  suggests  that  large  normal 
stress  may  be  transferred  into  a  fiber  even  if  the  interface  shear  stress  is  very  low. 
Therefore,  apart  from  fiber  breaks,  fiber  pullout  is  not  expected,  and  the  integrity  of  the 
fibrous  plies  in  the  axial  direction  should  not  be  impaired  in  the  presence  of  fiber  sliding. 

Guided  by  results  of  our  related  studies,  which  suggest  that  large  diameter  fibers 
release  a  substantial  amount  of  energy  during  debonding  and  thus  experience  nearly  total, 
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instantaneous  separation,  we  represented  the  interface  failure  under  normal  or  shear 
stresses  by  reducing  certain  fiber  stress  components  to  zero  through  abrupt  changes  in  the 
fiber  elastic  moduli.  In  this  part,  the  fiber  was  regarded  as  transversely  isotropic,  with 
moduli  equal  to  isotropic  material  magnitudes  (see  the  Appendix).  If  the  interface  normal 
stress  component  reached  the  selected  strength  magnitude  of  50  MPa,  then  the  transverse 

f  f 

shear  and  Young’s  moduli  GT  and  ET  were  reduced  to  zero.  Similarly,  if  the  interface 
longitudinal  shear  stress  exceeded  the  selected  80  MPa  shear  strength,  the  longitudinal 

f 

shear  modulus  GL  was  reduced  to  zero.  This  approach  has  been  implemented  in  the 
transformation  analysis  method.  The  ply  elastic  moduli  in  (22)  and  (23)  have  been  found 
by  the  Mori— Tanaka  method.  After  fiber  moduli  reductions,  the  ply  stiffnesses  and 
compliances,  as  well  as  the  concentration  factors,  were  reevaluated  and  introduced  in  the 
governing  equations  (21).  More  advanced  stages  of  damage,  such  as  matrix  cracking  and 
fiber  breaks,  were  not  considered  in  the  present  analysis. 

3.2.4  Finite  element  scheme 

We  now  consider  the  finite  element  solution  for  an  assemblage  of  M  elements 
interconnected  at  N  nodal  points  subjected  to  concentrated  loads  p^t),  i  =  1,  2,  ..  ,  N,  and 
uniform  temperature  variation  0m(t),  m  =  1,  2,  ..  ,  M.  Incorporating  the  transformation 
field  analysis  described  above  in  the  initial  strain  finite  element  formulation  given  by 
Zienkiewicz  and  Cormeau  [31],  Bahei— El— Din  [17]  developed  a  finite  element  scheme  for 
multiphase  materials  and  structures.  For  two— phase  materials  with  the  rate— dependent 
behavior  described  in  Section  2,  the  algorithm  for  computing  the  time— rate  of  the  nodal 
displacements,  i  =  1, 2,  ..  ,  N,  element  overall  strains  and  stresses,  em,  aw 
m  =  1,  2,  ..  ,  M,  and  phase  strains  and  stresses,  ev  av  r  =  1,  2,  is  outlined  for  infinitesimal 
strains  as  follows  [17]: 
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1.  Compute  the  overall  stiffness  matrix  K  at  current  temperature: 

Kij  =  X  /  tA^T  L»w  A"  dv» 

m=  1 

where  the  coefficient  matrix  Am  defines  the  element  average  strains,  cm,  in  terms  of 
the  nodal  displacements,  am;  cm  =  Am  am.  The  form  of  Am  depends  on  the  shape 
functions  selected  to  represent  the  displacement  field  in  the  elements.  If  the  strains 
in  element  m  is  required  to  be  uniform,  as  for  example  in  idealization  of  two— phase 
materials  by  averaging  models,  then  the  displacement  shape  functions  must  be 
linear  in  the  coordinates  Xj,  j  =  1,  2,  3,  and  matrix  Am  is  constant. 

2.  For  elements  m  =  1,  2,  ,  M: 

a)  From  the  constitutive  law  presented  in  eqs.  (1)— (3),  and  (7),  compute  the 
transformation  strain  rate  at  temperature  6m  for  phases  r  =  1,  2  as  the  sum  of 
viscoplastic  and  thermal  components; 


b) 


.(rm)  (rm)  .  < rm)  .  .  (rm)  "  f  ..  p(  en 

M ij  “  l  aij  +  (  ^  (Tkl  ]  6m  +  [k(0)  R  J 

Arm  =  {Mil  ^22  ^33  ^^23  ^31  2/i12}rm  • 

Compute  the  inelastic  strain  rate  at  time  t  for  element  m; 


rm 


V 


( rm) 

ij 


Atn  =  Ct  Aim  +  c2  Bl(eJ  km- 
c)  Compute  the  contribution  of  element  m  to  nodal  load  rates; 

f ?  =  -/  [a5t  uu  K, dv.  • 

V» 

where  Lm  is  the  overall  stiffness  matrix  of  element  m  evaluated  at  the  current 
temperature  of  the  element. 
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3.  Assemble  current  rates  of  nodal  loads: 


M 


m=  1 


4.  Compute  rates  of  nodal  deflections  at  time  t: 

a  =  K"1  f 

5.  For  elements  m  =  1,  2,  ,  M; 

a)  Compute  overall  strain  rate;  cm  =  Am  am , 

b)  Compute  overall  stress  rate;  am  =  Lm(0m)  (  em  -  jt*  ) , 

c)  For  phase  r  =  1,  2,  of  element  m; 

-  Compute  local  strain  rate;  er  =  Ar(0m)  <•  +  Drl(0m)  +  D t2(6J  ]i2  , 

—  Compute  local  stress  rate;  err  =  Lr(0m)  (  cr  —  At )  * 

The  above  rate  equations  are  assembled  into  a  system  of  first  order  differential 
equations  (ODE)  of  the  form 

yj  (0  =  gj  ( t,  yv  y2,  -  ,  yR )  • .  (27) 

The  unknown  functions  gj,  j  =  1,  2, ..  ,  R,  are  identified  with  the  nodal  displacements  ait 
i  =  1,  2,  ..  ,  N,  the  element  stress  <rm  and  strain  cm,  m  =  1,  2, ,  M,  the  phase  stress  ox  and 
strain  er,  r  =  1,  2,  and  any  internal  variables  required  to  define  the  eigenstrain  rate  f^,  e.g. 
the  scalar  function  R  as  in  eq.  (6).  Assuming  elastic  response  of  the  phases  in  the  initial 
state,  the  rate  equations  can  be  integrated  for  a  specified  time  period  using  an  ODE  solver 
for  stiff  differential  equations  which  are  normally  encountered  in  viscoplastic  response 
modeled  with  the  power  law  adopted  in  eq.  (7).  Implementation  of  the  finite  element 
algorithm  in  the  VISCOPAC  program  [18]  used  the  Gear  method  for  stiff  equations 
[32,  33]. 
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3.3  Interpretation  of  Isothermal  Static  and  Fatigue  Data  of  Sigma/Timetal  21S  Laminates 

Tk.  fatigue  experiments  described  in  Section  2  are  simulated  using  the 
analytical/numerical  approach  described  above.  The  computations  were  performed  for 
several  specimens  of  the  Sigma/Timetal  21S,  (0)4,  (0/90)s,  and  (0/±45/90)s  laminates 
tested  in  the  experimental  program.  Fabrication  by  hot  isostatic  pressing,  which  has  a 
significant  effect  on  the  overall  response,  was  modeled  first.  Subsequent  loading  consisted 
of  axial  tension/tension  load  cycles  at  a  specified  amplitude  and  frequency,  after  a  uniform 
change  of  temperature  from  21°C  to  650°C  in  6  minutes.  Assuming  that  the  overall  stress 
and  strain  were  uniform  in  the  gauge  section  the  specimen,  where  the  strain  measurements 
were  made,  the  response  of  the  laminates  under  the  applied  load  was  determined  from  the 
classical  laminated  plate  theory  for  in— plane  loading,  modeled  with  the  finite  element 
method. 

A  schematic  of  the  finite  element  model  for  laminates  is  shown  in  Fig.  26  for  the 
(0/±45/90)s  layup.  Each  unidirectional  ply  is  represented  by  one  hexahedral  element.  The 
overall  elastic  moduli  of  the  plies,  and  the  phase  mechanical  and  transformation 
concentration  factors  were  derived  from  the  Mori— Tanaka  model  of  the  fibrous 
microstructure.  As  described  in  Section  3.3,  interface  failure  by  fiber  separation  was 
assumed  to  take  place  when  the  maximum  interface  normal  stress  component  reached  a 
selected  normal  bond  strength  magnitude  <rb,  and  by  fiber  sliding  when  the  maximum 
interface  longitudinal  shear  stress  exceeded  a  selected  shear  bond  strength  rb.  In  analysis 
of  the  Sigma/Timetal  21S  system,  <rb  =  50  MPa  and  rb  =  80  MPa  were  assumed  for  the 
interface  strength  magnitudes  against  separation  and  sliding,  respectively. 

3.3.1  Hot  isostatic  pressing 

Fabrication  of  the  Sigma/Timetal  21S  laminates  by  hot  isostatic  pressing  was 
conducted  at  899°C  (1650°F)  and  103.5  MPa  (15  ksi)  hydrostatic  pressure  for  2  hours.  In 
the  analytical  simulation,  we  assumed  that  the  laminate  is  free  of  internal  stresses  at  the 
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fabrication  temperature  where  the  matrix  yield  stress  is  very  low,  Table  6.  Subsequent 
cooling  to  room  temperature  was  completed  in  3  minutes  and  followed  by  aging  at  621°C 
(1150°F)  for  8  hours.  Figure  27  shows  the  temperature  and  load  history  applied  during 
fabrication  and  aging.  We  note  that  the  actual  aging  period  of  8  hours,  during  which  the 
titanium  matrix  remained  elastic,  has  been  reduced  to  2  hours  in  Fig.  27. 

The  average  residual  stresses,  referred  to  the  local  coordinate  system  Xj,  j  =  1,  2,  3, 
Fig.  26,  caused  by  fabrication  in  the  fiber  and  matrix  of  the  (0)4  unidirectional  laminate 
are  shown  in  Table  7  after  cooldown  to  room  temperature  and  subsequent  reheating  to  the 
test  temperature  (650°C).  The  fiber  supports  a  compressive  axial  residual  stress  aw  which 
increases  the  apparent  strength  of  the  composite  under  subsequent  axial  tensile  loading. 
The  fiber  also  supports  a  compressive  transverse  stress  a22  after  fabrication.  A  more 
detailed  study  of  the  local  fields  computed  in  a  composite  cylinder  model  of  a  unidirectional 
titanium  matrix  composite  [34]  indicated  that  the  radial  residual  stress  at  the  fiber 
interface  was  indeed  compressive,  and  thus  supportive  of  interface  bonding.  This  suggests 
that  tensile  hoop  stress  is  present  in  the  matrix  at  the  fiber  interface.  This  may  promote 
radial  cracking  in  the  matrix  and  open  the  way  for  environmental  degradation  of  the  fiber, 
particularly  at  high  temperature.  The  magnitude  of  the  residual  stresses  may  be  adjusted 
somewhat  by  varying  the  cooling  rate  and  the  magnitude  of  the  HIP  pressure  [34,  35]. 

The  local  stresses  caused  in  the  plies  of  quasi— isotropic  (0/90)s  and  (0/±45/90)s 
laminates  by  the  thermomechanical  loads  applied  in  the  HIP  process  are  expected  to  be 
independent  of  the  fiber  orientation  [36,  37].  Table  8  shows  the  residual  stresses  found  in  a 
typical  ply  of  these  laminates  and  in  the  fiber  and  matrix  phases  after  cooldown  to  room 
temperature  and  subsequent  reheating  to  the  test  temperature  (650°C).  As  in  the 
unidirectional  composite,  the  residual  axial  fiber  stress  an  is  compressive  and  thus  helps  to 
reduce  the  tensile  stress  magnitude  under  subsequent  tensile  mechanical  loading.  However, 
a  tensile  residual  stress  <r22  =  15  MPa  is  present  in  the  fiber  at  650°C,  which  may  promote 
debonding  under  subsequent  axial  tension  loading.  Of  course,  under  the  selected  value  of 
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the  interface  normal  bond  strength  <rb  =  50  MPa,  the  fiber  interface  remains  intact  at  the 
test  temperature  of  650°C. 

The  overall  axial  strains  were  measured  experimentally  at  650°C  as  0.40,  0.46,  and 
0.45%  in  the  (0)4,  (0/90)s,  and  (0/±45/90)s  laminates,  respectively,  Table  1.  They 
compare  well  with  the  computed  strains  of  0.377%  for  the  (0)4  unidirectional  composite, 
and  0.424%  for  each  of  the  quasi— isotropic  (0/90)s  and  (0/±45/90)s  laminates. 

3.3.2  Static  loading 

In  the  actual  fatigue  experiments,  Section  2.4,  mechanical  loading  to  the  mean 
stress  was  applied  in  one  minutes  and  kept  constant  for  one  minute  before  application  of 
the  load  cycles.  Consequently,  the  mean  stress  load  rates  varied  according  to  the  stress 
range  selected  for  the  various  layups.  In  the  numerical  simulation  of  this  load  step,  fiber 
debonding  and  sliding  took  place  in  the  off-axis  plies  under  overall  axial  tensile  stress 
magnitudes  lower  than  the  mean  stress,  with  interface  normal  strength  <rb  =  50  MPa  and 
shear  strength  rb  =  80  MPa.  Figures  28  and  29  show  variation  of  the  transverse  normal 
stress  and  longitudinal  shear  stress  in  the  fully  bonded  fiber  of  the  off-axis  plies  of  the 
Sigma/Timetal  21S  (0/90)s  and  (0/±45/90)s  laminates  as  a  function  of  the  applied  tensile 
axial  stress  at  the  two  applied  loading  rates.  The  offset  stress  found  at  zero  overall  load 
represents  the  residual  stress  after  fabrication,  Table  8.  At  the  selected  <rb  =  50  MPa,  fiber 
debonding  initiates  under  an  overall  axial  stress  of  about  25  MPa.  Similarly,  initiation  of 
interface  sliding  in  the  45°  plies  at  rb  =  80  MPa  takes  place  while  the  overall  axial  stress  is 
less  than  100  MPa.  However,  at  these  overall  load  magnitudes,  the  rate  of  loading  has  no 
significant  effect  on  the  local  stresses. 

To  illustrate  the  extent  of  interface  damage  which  is  likely  to  take  place  in  the 
Sigma/Timetal  21S  specimens  during  the  fatigue  tests,  we  computed  the  local  stresses  in  a 
(0/±45/90)s  laminate  loaded  by  an  overall  axial  stress  at  2.75  MPa/s.  This  rate 
corresponds  to  loading  to  the  mean  stress  of  the  300—30  MPa  stress  range  in  one  minute. 
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able  9  show,  ,he  competed  local  s.tesses  in  the  plies  and  in  the  fiber  and  matrix  phases  at 
the  onset  of  mterface  failure  in  the  individual  plies,  together  with  the  corresponding  overall 
stress  magnitudes.  Interface  failure  initiates  in  the  90^egree  fiber  when  the  transverse 
tensile  stress  reaches  the  selected  interface  strength  of  50  MPa,  at  an  overall  axial  stress  of 
only  18  MPa.  Subsequent  loading  to  26  MPa  caused  fiber  debonding  in  the  degree 
Plies,  followed  by  debonding  of  the  0-degree  fiber  at  59  MPa.  Finally,  the  magnitude  of 

longitudinal  shear  stress  in  the  fiber  of  the  degree  plie,  reached  the  selected  shear 
strength  of  80  MPa  at  overall  axial  stress  of  65  MPa. 

Under  subsequent  loading,  the  debonded  fibers  in  the  off-axis  plies  as  well  as  the 
0-degree  ply  fibers,  support  only  axial  normal  stresses.  Consequently,  the  stiffness  of  the 
off-axis  pl.es  in  the  loading  direction  is  reduced  substantially,  and  the  axial  stress 
concentration  factor,  given  by  the  ratio  of  the  local  to  overall  axial  stress  increments 
increases  in  the  Degree  layer.  For  example,  as  seen  from  Table  9,  this  ratio  increased 
from  2.8  a.  the  overall  stress  of  18  MPa,  where  all  fibers  are  still  perfectly  bonded,  4  2  at 
66  MPa,  where  the  fiber  were  debonded  in  all  plies.  The  corresponding  ratios  for  the  fiber 
stress  in  the  0-degree  p.y  are  5.4  and  9.4,  which  indicate  a  substantial  increase  in  the  fiber 
stress  in  the  axial  direction  due  to  interface  damage  of  the  off-axis  plies.  At  165  MPa,  the 
overall  axial  strain  was  computed  at  0.72%,  sfightly  above  the  measured  value  of  0.64%’. 

3.3.3  Isothfrrmal  load  cycles 

The  Local  stresses  computed  for  the  Sigma/Timetal  21S  (0/*45/90)s  laminate 
subjected  to  axial  tension/tension  load  cycles  at  650°C,  with  maximum  stress 

S“"  ~  300  MPa'  Smi”/S»“  =  01>  lnd  a  P™»d  of  lOOOs/cycle  are  summarized  in 
Figs.  30-39.  Figure  30  shows  the  load  history  during  the  eleven  cycles  sustained  by  the 

spectmen  in  the  actual  experiments,  Section  2.4.  Figure  3  shows  the  average  stresses 
generated  in  the  O’  ply,  and  Figs.  32  and  33  provide  the  local  stress  averages  in  the  fiber 
and  matrix.  The  ply  and  local  stresses  in  the  *45-degree  plies  are  shown  in  Figs.  34-36, 
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and  the  stresses  in  the  90— degree  ply  are  presented  in  Figs.  37—39.  The  ratcheting 
behavior  found  iii  the  stress  histories  is  caused  by  the  viscoplastic  response  of  the  titanium 
matrix,  modeled  by  the  constitutive  equations  given  in  Section  3.1.  In  general,  the  matrix 
stresses  are  relaxed  and  the  fiber  axial  stress  in  all  plies  is  enhanced  under  the  applied 
tension/tension  load  cycles.  Under  slow  load  cycles  applied  at  0.001  Hz,  a  nearly  steady 
state  response  is  reached  after  only  few  cycles. 

The  cyclic  stress  maxima  and  minima  computed  in  the  first  cycle  and  in  the  failure 
cycle  11  are  shown  in  Table  10  under  the  cyclic  loading  conditions  indicated  in  Fig.  30  for 
the  slow  fatigue  cycles  at  0.001  Hz  (Figs.  31-39).  The  stresses  computed  for  a  fast  cycle  at 
0.1  Hz  under  the  same  loading  conditions  are  also  indicated.  Another  test  sample  subjected 
to  the  same  load  amplitudes,  but  at  0.1  Hz,  survived  207  cycles,  Table  2.  In  this  case,  the 
computed  local  stresses  reached  a  steady  state  after  90  cycles.  The  local  stresses  were  then 
extrapolated  to  the  failure  cycle  using  a  logarithmic  function  of  the  cycle  number; 
<r  =  A  +  B  log  N,  where  N  is  the  cycle  number,  and  A  and  B  are  fitting  parameters  to  be 
computed  from  the  data  at  the  steady  state.  For  example,  using  the  computed  data  at 
cycles  91,  92,  and  93,  at  which  the  local  stresses  stabilize,  the  average  axial  stress  <Jn  in  the 
0°  ply  is  extrapolated  to  the  failure  cycle  207  using  the  above  logarithmic  equation  with 
A  =  735.7,  and  B  =  63.2. 

The  rate  effect  on  the  local  stresses  is  evident  in  the  results  given  in  Table  10  for  the 
first  cycle.  For  example,  the  matrix  axial  normal  stress  <ju  =119  MPa  when  computed  in 
the  0— degree  ply  at  Smax=  300  MPa  under  0.001Hz.  At  0.1  Hz  and  Smax=  300  MPa, 
cr^  =  206  MPa.  In  contrast,  the  axial  stress  in  the  0°  fiber  is  elevated  at  the  slow  loading 
rate.  The  maximum  fiber  stress  found  at  failure,  however,  appears  to  be  independent  of 
the  loading  frequency.  Specifically,  the  axial  stress  computed  in  the  0°  fiber  at  failure  is 
1981  MPa  under  0.001  Hz,  and  1988  MPa  under  0.1  Hz.  This  suggests  that  failure  of  the 
laminate  is  controlled  by  the  strength  of  the  fiber  in  the  0°  ply,  and  that  the  fatigue  life  is 
dependent  on  the  rate  at  which  the  fiber  tensile  strength  is  reached. 
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Both  the  average  transverse  and  longitudinal  shear  stresses  in  the  0°  ply  matrix  are 
low  in  plies  with  debonded  fibers,  and  decrease  with  cycling.  Also,  the  maximum  axial 
stress  in  the  ±45°  fiber  an  =  687  MPa,  and  remains  low  in  the  ±45°  fiber. 

A  failure  Criterion  may  be  postulated  for  laminates  subjected  to  cyclic  axial  loading 
based  on  the  preceding  discussion  of  local  stresses.  In  particular,  a  one— parameter  criterion 
based  on  the  tensile  strength,  or  ultimate  strain  of  the  0°  ply,  or  its  fiber  may  be 
representative  of  the  observed  fatigue  behavior.  The  overall  axial  strains  measured  at 
failure  for  the  Sigma/Timetal  21S  laminates  tested  in  the  experimental  program  under 
tension/tension  load  cycles  applied  at  650°C  are  shown  in  Table  2. 

Although  the  measured  fatigue  life  of  the  tested  laminates  varied  with  the  load 
frequency,  Fig.  23,  the  overall  strains  measured  at  failure  were  not  frequency— dependent. 
Moreover,  the  measured  failure  strains  appear  to  be  independent  of  both  the  layup  and  the 
stress  amplitude.  This  is  shown  in  Figs.  40  and  41  where  the  overall  strain  at  failure  is 
plotted  as  a  function  of  the  number  of  cycles.  The  measured  strains,  represented  by  the 
closed  symbols,  fall  within  a  narrow  range  in  the  order  of  1.0  ±  0.1  %  for  both  the  fast  and 
the  slow  load  cycles. 

Predictions  of  the  overall  strains  are  indicated  by  the  open  symbols  in  Figs.  40 
and  41.  The  computed  strains  are  consistent  with  the  measured  values  and  fall  in  the 
experimental  strain  range.  This  suggests  that  the  fatigue  life  of  composite  laminates  can  be 
estimated  based  on  the  computed  total  overall  strain.  For  a  ductile  matrix  with 
rate— dependent  behavior,  such  as  the  titanium  matrix  examined  here,  the  number  of  cycles 
for  reaching  the  critical  strain  magnitude  is  a  function  of  the  load  frequency. 
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6.  CONCLUSIONS 


The  experimental  results  show  that  at  650°C,  viscoplastic  deformation  of  the 
Timetal  21S  matrix  can  be  detected  even  at  relatively  low  applied  stresses.  At  both  650°C 
and  21°C,  reduction  of  unloading  elastic  moduli  of  the  (0/90)s  and  (0/±45/90)s  laminates 
was  observed  after  either  cyclic  or  static  loading  to  relatively  low  tensile  stresses.  This 
indicates  onset  of  damage.  In  static  tests  at  room  temperature,  interface  decohesion  and 
fiber  sliding  appear  to  be  the  principal  damage  modes.  In  fatigue  tests  at  650°C,  matrix 
cracking  is  also  present.  The  cyclic  loading  rate  has  no  effect  on  endurance  of  the  (0/90)s 
and  (0/±45/90)s  laminates  loaded  by  constant  amplitude  tension  at  650°C.  In  the  (0)4 
laminate,  higher  endurance  limits  were  observed  at  0.1Hz  that  at  0.001  Hz.  However, 
regardless  of  rate  and  layup,  the  total  strain  at  failure  under  both  static  and  cyclic  loading 
at  650°C,  was  measured  at  1±0.1%.  In  room  temperature  static  tests,  the  failure  strain  was 
measured  at  0.5-0. 6%, 

Micromechanical  modeling  of  viscoplastic  composite  laminates  with  the 
transformation  field  analysis  (TFA)  and  the  finite  element  method  has  been  described  and 
applied  to  Sigma/Timetal  21S  laminates  using  the  Mori— Tanaka  model  for  unidirectional 
composites.  The  computed  results  provide  a  qualitative  evaluation  of  the  local  stresses 
under  fabrication  and  subsequent  loading,  and  point  to  possible  causes  of  fatigue  failure  in 
(0/90)s  and  (0/±45/90)s  laminates. 

Under  constant  temperature,  the  fiber  stresses  in  the  0°  plies  increase  with  the 
number  of  cycles,  until  a  certain  overloading  situation  is  reached.  The  rate  of  change  of 
the  fiber  stress  with  cycles  is  reflected  by  the  ratcheting  rate  of  the  viscoplastic  titanium 
matrix,  particularly  at  high  temperature,  and  is  a  function  of  the  loading  frequency.  The 
axial  stresses  estimated  in  the  0°  fiber  under  identical  isothermal  loading  conditions,  but 
different  frequencies  are  equal.  Substantial  interface  damage  is  expected  to  take  place  in 
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the  fiber  of  the  off-axis  plies,  but  the  matrix  stresses  are  small  and  may  not  cause 
significant  matrix  cracking. 

The  fatigue  data  measured  for  the  Sigma/Timetal  21S  laminates  under  various 
stress  amplitudes,  and  two  frequencies  which  differ  by  two  orders  of  magnitude,  suggest  a 
one— parameter  failure  criterion  based  on  a  1  ±  0.1  %  ultimate  strain  of  the  0°  ply,  or  its 
fiber,  at  650°C.  This  maximum  measured  strain  is  closely  approximated  by  the  strain 
computed  at  the  failure  cycle  using  the  micromechanical  TFA/FE  modeling  approach. 
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APPENDIX 


Suppose  that  the  fiber  and  matrix  moduli  and  coefficients  of  thermal  expansion  are 
represented  by  continuous  functions  of  temperature  6,  derived  from  experimental  data. 
After  a  uniform  change  AO,  from  a  reference  temperature  60,  applied  under  a  uniform  stress 
aT,  an  unconstrained  phase  deforms  uniformly, 

0 

er  =  Mr(0)  aT  +  J  m.(0)d0.  (A-l) 

K 


If  the  temperature  changes  under  an  applied  uniform  strain,  the  stress  is: 


*r  =  Lt(0) 


/  “,(0)  » 


(A— 2) 


where  r  =  f,  m  denotes  the  phase,  Lr(  6)  and  Mr(<?)  =  [Lr(0)]'1  are  temperature  dependent, 
(6x6)  phase  stiffness  and  compliance  matrices,  m^tf)  is  the  (6x1)  thermal  strain  vector  of 
expansion  coefficients,  their  particular  forms  depend  on  the  elastic  symmetry  of  the  phase. 

For  a  transversely  isotropic  solid,  with  xt  as  the  axis  of  rotational  symmetry,  and 
compliance  matrix  M,  the  constitutive  relation  e  =  M  a  is  written  using  the  engineering 
moduli 
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The  stiffness  L  =  M'1  is  best  written  in  terms  of  Hill’s  [38]  moduli, 

k  =  -  [1/G23  — 4/E22  +  4^,/Ej'1 ,  1  =  2k  vu 

n  =  Eu  +  4k  u2n  =  Ea+  P/k  ,  m  =  G23 , 


(A— 3) 


P  —  G12  • 


The  constitutive  relation  a=L  e  then  is 
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Estimates  of  the  overall  moduli  L  of  heterogeneous  solids  by  the  Mori— Tanaka 
method  are  given  in  references  [22,  23],  Also,  the  following  formulas  for  the  overall  Hill’s 
moduli  of  a  fibrous  composite  in  terms  of  the  local  moduli  and  the  fiber  volume  fraction  cf 
were  derived  by  Chen  et  al.  [39]  for  the  Mori-Tanaka  method; 


2cf  Pm  Pf  +  cm(pm  pf  +  pj) 

P  =  .  .  > 

2cf  Pm  +  cm(Pf  +  Pm) 

(A— 5) 

mrn  mf  (km  +  2mrn)  +  km  mm  ( c  f  “f  +  cm  mrn) 
m  =  , 

km  mrn  +  (k m  +  2mrn)  (  cf  mrn  +  cm  mf) 

(A— 6) 

cf  kf  (km  +  mj  +  cm  km  (kf  +  rnm) 

k  = - , 

cf  (km  +  mm)  +  cm  (kf  +  rnm) 

(A-7) 

cf  l{  (km  +  rnm)  +  cm  4  (kf  +  rnm) 

1= - , 

Cf  (km  +  mrn)  +  cm  (kf  +  mrn) 

(A-8) 

—  ^m 

n  =  cm  nm  +  cf  nf  +  {1  cf  If  cm  4) 

kf  -  km 

(A— 9) 

42 


I 

I 


Table  1  Results  of  static  tests  performed  on  Sigma/Timetal  2 IS  laminates  at  650°C. 


Thermal 

Stress 

Elastic 

Ultimate 

Ultimate 

Layup 

Strain 

Rate 

Modulus 

Strength 

Strain 

(%) 

(MPa/s) 

(GPa) 

(MPa) 

(%) 

(0)4 

0.40 

0.278 

1.205 

27.80 

157 

819a 

1.040 

(0/90)s 

0.46 

0.278 

100 

408a 

1.076 

27.80 

125 

482a 

0.946 

(0/±45/90)s 

0.45 

114 

348b 

1.227 

125 

431a 

0.962 

a  Plate  A.  b  Plate  B. 


Table  2  Experimental  data  for  isothermal  fatigue  tests  of  Sigma/Timetal  21S 
laminates  (R  =  0.1,  0  =  650°C). 


Laminate 

f 

=  0.001 

Hz 

f 

=  0.1  Hz 

Smax  (MPa) 

Nf 

fult  (%) 

Smax  (MPa) 

Nf 

eult  (%) 

700ba 

1018* 

0.936 

850b 

1223 

1.165 

600b(3 

1750* 

0.960 

775b 

558 

1.079 

575b 

17 

0.939 

750a 

5700 

0.993 

562a 

2 

0.969 

725b 

2604 

1.031 

500a 

17 

1.074 

700b 

10191 

1.021 

450a 

208 

1.110 

425a 

2020 

1.118 

(0/90)s 

325a 

14 

0.989 

350a 

240 

0.960 

300a 

1182* 

0.904* 

325a 

15906 

0.984 

300a 

113750* 

0.884* 

200a 

110000* 

0.755* 

(0/±45/90)s 

300a 

11 

1.150 

300b 

207 

1.110 

250a 

215 

1.170 

250b 

2850 

1.070 

250b 

1201* 

1.152* 

250b 

12800 

1.172 

200b 

105515* 

0.941* 

200b 

18313* 

0.960* 

a  Plate  A.  b  Plate  B. 

«  527  cycles  at  600-*60  MPa  +  235  at  650-^65  +  251  at  675^67.5  +  5  at  700^70. 
*  Survived  the  number  of  cycles  and  strain  indicated. 
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Table  3  Unloading  elastic  modulus  (GPa)  measured  in  qusai-static  loading/unloading 
tests  performed  on  (0/90)s  Sigma/Timetal  21S  laminates  at  room  temperature. 


St  r ess 
(MPa) 

Reference  Specimen 

Specimen  —  1 

Specimen  —  2 

100 

160 

150 

159 

200 

150 

147 

150 

300 

144 

147 

350 

141 

400 

136 

450 

135 

Isothermal  fatigue  at  650°C: 

475 

134 

500 

134 

200  -*-20  MPa 

300  -4  30  MPa 

550 

132 

0.1  Hz 

0.1  Hz 

600 

130 

110000  cycles 

113750  cycles 

650 

129 

Suit  =  650  MPa 

294  MPa 

350  MPa 

Table  4  Unloading  elastic  modulus  (GPa)  measured  in  qusai-static  loading/unloading 
tests  performed  on  (0/±45/90)s  Sigma/Timetal  21S  laminates  at  room  temperature. 


St  ress 
(MPa) 

Reference  Specimen 

Specimen  - 1 

Specimen  —  2 

100 

155 

139 

140 

200 

150 

135 

135 

300 

134 

125 

131 

350 

125 

119 

131 

400 

119 

127 

Isothermal  fatigue 

at  650°C: 

250  -»  25  MPa 

200  -*  20  MPa 

0.001  Hz 

0.1  Hz 

1201  cycles 

105000  cycles 

Suit  =  437  MPa 

400  MPa 

400  MPa 
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Table  5  Elastic  moduli  and  coefficients  of  thermal  expansion 


Sigma  fiber 

Timetal  21S  matrix 

8 

(»C) 

E 

(GPa) 

V 

a 

(io-6/°c) 

8 

(°C) 

E 

(GPa) 

V 

a 

(10'6/°C) 

21 

393 

0.25 

3.56 

21 

116 

0.338 

8.62 

93 

390 

0.25 

3.56 

150 

116 

0.347 

8.72 

204 

386 

0.25 

3.62 

316 

104 

0.358 

9.23 

426 

378 

0.25 

3.91 

427 

103 

0.363 

9.87 

538 

374 

0.25 

4.07 

482 

96 

0.365 

10.30 

649 

370 

0.25 

4.27 

565 

86 

0.365 

11.10 

760 

365 

0.25 

4.41 

621 

80 

0.365 

11.50 

871 

361 

0.25 

4.57 

650 

64 

0.365 

11.80 

704 

53 

0.374 

12.00 

760 

50 

0.382 

12.50 

Table  6  Viscoplastic  parameters  for  Timetal  21S 


8 

(«C) 

Y 

(MPa) 

¥ 

(MPa) 

Ho 

(MPa) 

h 

(MPa) 

m 

k 

(MPa) -p/s 

P 

21 

925 

1000 

72 

4024 

1.365 

2.10E— 21 

8.59 

150 

925 

1000 

72 

4024 

1.365 

2.10E-21 

8.59 

316 

645 

785 

146 

5432 

1.355 

2.92E-16 

6.52 

427 

630 

770 

143 

5432 

1.365 

1.66E-14 

5.62 

482 

505 

615 

103 

3165 

1.478 

5.38E-13 

6.35 

565 

190 

365 

256 

6311 

1.211 

6.07E-12 

5.52 

621 

85 

185 

149 

2618 

1.200 

6.05E-10 

3.58 

650 

30 

87 

42 

1043 

1.000 

1.76E— 08 

2.00 

704 

24 

53 

45 

829 

1.000 

1.50E— 08 

2.00 

760 

9 

20 

15 

250 

1.000 

1.00E— 08 

2.00 

Table  7  Nonzero  local  stress  averages  (MPa)  caused  in  the  fiber  and  matrix  by 
hot  isostatic  pressing  of  a  unidirectional  Sigma/Timetal  21S  composite 


Temperature  Fiber  Stresses  Matrix  Stresses 


(»C) 

0n 

cr 22 

033 

0u 

022 

033 

21 

-770 

-228 

-228 

553 

158 

158 

650 

-136 

-50 

-50 

94 

35 

35 

Table  8  Nonzero  local  stress  averages  (MPa)  caused  by  hot  isostatic  pressing  in  the  plies  of 
Sigma/Timetal  21S  quasi-isotropic  (0/90)s  and  (0/±45/90)s  laminates 


Temperature 

Ply  Stresses 

Fiber  Stresses 

Matrix  Stresses 

(«C) 

011 

022 

0n 

022 

033 

0n 

022  033 

21 

-157 

157 

-1095 

-52  - 

-232 

492 

302  161 

650 

-55 

55 

-293 

15 

-53 

110 

82  36 
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Table  9  Local  stress  averages  (MPa)  in  plies  of  a  Sigma/Timetal  21S  (0/±45/90)s  laminate 

caused  by  overall  static  axial  loading  at  650°C. 


Ply 

Locala 

Under  Overall  Axial  Stress  (MPa) 

Stress 

Component 

18 

90°  ply 
debonding 

26 

±45°  plies 
debonding 

59 

0°  ply 
debonding 

65 

±45°  plies 
sliding 

0° 

<Uip 

-5 

39 

210 

217 

<?22  p 

54 

59 

85 

40 

<Uif 

-195 

-109 

256 

320 

(T22f 

15 

20 

50 

0 

<7'33f 

-53 

-55 

-65 

0 

<Ulm 

127 

141 

177 

145 

022m 

82 

86 

110 

67 

033m 

37 

38 

45 

0 

45“ 

<ruP 

-33 

-10 

53 

77 

022P 

69 

84 

47 

50 

0"12P 

15 

26 

54 

54 

an{ 

-255 

-214 

-30 

21 

022* 

32 

50 

0 

0 

<7  33f 

-54 

-56 

0 

0 

<U2f 

20 

34 

78 

80 

<Uim 

121 

131 

110 

115 

<T22m 

95 

108 

80 

85 

^33m 

37 

39 

0 

0 

ai2m 

11 

19 

37 

37 

CO 

o 

o 

<7\ lp 

-61 

-82 

-77 

-58 

022p 

84 

43 

61 

65 

ffnf 

-314 

-325 

-319 

-281 

<^22f 

50 

0 

0 

0 

(T33f 

-54 

0 

0 

0 

an* 

115 

85 

90 

97 

<?22m 

107 

72 

103 

110 

&23m 

38 

0 

0 

0 

a  o-.jp  =  ply  stress,  <7ijf  =  fiber  stress,  crijm  =  matrix  stress;  xi  =  axial  direction 
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Table  10  Local  stress  averages  (MPa)  in  the  plies  of  a  Sigma/Timetal  21S 
(0/±45/90)s  laminate,  caused  by  axial  load  cycles 
(Sraax  =  300  MPa,  R  =  0.1,  9  =  650°C). 


Ply 

Plv 

f  = 

0.001  Hz 

f  = 

0.1  Hz 

riy 

Q+rooo 

cycle  1 

cycle  lla 

cycle  1 

cycle  207a 

O II  Coo 

max. 

min. 

max. 

min. 

max. 

min. 

max. 

min. 

0” 

<TU  P 

786 

121 

840 

165 

732 

86 

882 

231 

<T22p 

23 

7 

-11 

-18 

31 

27 

14 

10 

anf 

1749 

333 

1981 

513 

1492 

243 

1988 

713 

anm 

119 

-25 

51 

-77 

206 

-22 

117 

-103 

<?22m 

40 

13 

-19 

-31 

53 

47 

23 

18 

45° 

crnP 

257 

29 

302 

65 

239 

8 

270 

38 

CT22P 

43 

0 

1 

-33 

62 

24 

15 

-24 

tfl2p 

24 

-11 

18 

-16 

36 

-10 

4 

-41 

anf 

506 

65 

687 

208 

389 

-32 

613 

185 

(Tlim 

84 

3 

35 

-34 

135 

35 

34 

-63 

&  22m 

73 

1 

2 

-56 

106 

40 

23 

^1 

al2m 

40 

-18 

31 

-26 

60 

-17 

7 

-69 

CO 

o 

o 

0up 

-276 

-58 

-256 

-44 

-261 

-78 

-289 

-106 

0‘22P 

67 

-8 

21 

^5 

95 

22 

17 

-56 

an{ 

-739 

-205 

-609 

-100 

-717 

-309 

-756 

-343 

anm 

44 

44 

-11 

-5 

54 

81 

32 

59 

022m 

113 

-14 

35 

-76 

162 

37 

28 

-94 

a  Number  of  cycles  at  failure  found  experimentally  (Table  2) 
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FIGURE  CAPTIONS 


Figure  1.  Temperature  profile  obtained  for  a  Sigma/Timetal  21S  specimen 

in  a  succeptor  furnace 

Figure  2.  Applied  stress,  and  measured  total  strain  histories  in  a  multistep 

creep  test  for  a  Timetal  21 S  specimen  at  650°C 

Figure  3.  Creep  strains  computed  from  the  experimental  data  of  Fig.  2 

Figure  4.  Cyclic  loading  applied  to  a  Timetal  21S  specimen  at  650°C 

Figure  5.  Maximum  and  minimum  strains  measured  under  the  cyclic  loading  of  Fig.  4 

Figure  6.  Maximum  and  minimum  strains  produced  under  each  of  the  cyclic  stress 

ranges  of  Fig.  4 

Figure  7.  Tensile  stress-strain  data  measured  for  Sigma/Timetal  21S  laminates 

at  650°C  under  two  stress  rates 

Figure  8.  Total  axial  strain  histories  measured  in  multistep  creep  tests  for  two 

(0)4  Sigma/Timetal  21S  laminates  at  650°C 

Figure  9.  Total  axial  strain  history  measured  in  a  multistep  creep  test  for  a  (0/90)s 

Sigma/Timetal  21S  laminate  at  650°C 

Figure  10.  Total  axial  strain  history  measured  in  a  multistep  creep  test  for 
a  (0/±45/90)s  Sigma/Timetal  21S  laminate  at  650°C 

Figure  11.  Total  axial  strain  history  measured  for  a  (0)4  Sigma/Timetal  21S  laminate 
at  650°C  under  the  slow  load  cycles  shown  in  inset 

Figure  12.  Total  axial  strain  history  measured  for  a  (0)4  Sigma/Timetal  21S  laminate 
at  650°C  under  the  fast  load  cycles  shown  in  inset 

Figure  13.  Stress— axial  strain  response  measured  for  a  (0)4  Sigma/Timetal  21S 
laminate  under  the  slow  load  cycles  shown  in  inset  of  Fig.  11 

Figure  14.  Stress-axial  strain  response  measured  for  a  (0)4  Sigma/Timetal  21S 
laminate  under  the  fast  load  cycles  shown  in  inset  of  Fig.  12 

Figure  15.  Total  axial  strain  history  measured  for  a  (0/90)s  Sigma/Timetal  21S 
laminate  at  650°C  under  the  slow  load  cycles  shown  in  inset 

Figure  16.  Total  axial  strain  history  measured  for  a  (0/90)s  Sigma/Timetal  21S 
laminate  at  650°C  under  the  fast  load  cycles  shown  in  inset 

Figure  17.  Stress— axial  strain  response  measured  for  a  (0/90)s  Sigma/Timetal  21S 
laminate  under  the  slow  load  cycles  shown  in  inset  of  Fig.  15 

Figure  18.  Stress-axial  strain  response  measured  for  a  (0/90)s  Sigma/Timetal  21S 
laminate  under  the  fast  load  cycles  shown  in  inset  of  Fig.  16 
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Figure  19.  Total  axial  strain  history  measured  for  a  (0/±45/90)s  Sigma/Timetal  21S 
laminate  at  650°C  under  the  slow  load  cycles  shown  in  inset 

Figure  20.  Total  axial  strain  history  measured  for  a  (0/±45/90)s  Sigma/Timetal  21S 
laminate  at  650°C  under  the  fast  load  cycles  shown  in  inset 

Figure  21.  Stress— axial  strain  response  measured  for  a  (0/±45/90)s  Sigma/Timetal  21S 
laminate  under  the  slow  load  cycles  shown  in  inset  of  Fig.  19 

Figure  22.  Stress— axial  strain  response  measured  for  a  (0/±45/90)s  Sigma/Timetal  21S 
laminate  under  the  fast  load  cycles  shown  in  inset  of  Fig.  20 

Figure  23.  Maximum  stress— cycles  to  failure  (S— N)  data  measured  for 

Sigma/Timetal  21S  laminates  at  650°C  under  two  frequencies 

Figure  24.  Variation  of  secant  modulus  measured  in  two  (0)4  Sigma/Timetal  21S 

laminates  subjected  to  cyclic  loading  at  650°C 

Figure  25.  Maximum  strains  measured  at  failure  in  fatigue  tests  of  Sigma/Timetal  21S 

laminates  at  650°C 

Figure  26  Finite  element  model  of  a  (0/±45/90)s  laminate 

Figure  27  Thermomechanical  loading  applied  to  Sigma/Timetal  21S  laminates 

in  fabrication  by  hot  isostatic  pressing  and  aging 

Figure  28  Variation  of  the  transverse  normal  stress  in  fully  bonded  fibers  of 

a  Sigma/Timetal  21S  (0/90)s  laminate  with  the  overall  axial  stress 

Figure  29  Variation  of  the  transverse  normal  stress,  and  longitudinal  shear  stress  in  fully 
bonded  fibers  of  a  Sigma/Timetal  21S  (0/±45/90)s  laminate  with  the  overall  axial  stress 

Figure  30  Isothermal  cyclic  loading  applied  to  a  Sigma/Timetal  21S  (0/±45/90)s  laminate 

Figure  31  Stress  average  computed  in  the  0°  ply  of  a  Sigma/Timetal  21S  (0/±45/90)s 
laminate  subjected  to  HIP  and  isothermal  cyclic  loading 

Figure  32  Fiber  stress  average  computed  in  the  0°  ply  of  a  Sigma/Timetal  21S 
(0/±45/90)s  laminate  subjected  to  HIP  and  isothermal  cyclic  loading 

Figure  33  Matrix  stress  average  computed  in  the  0°  ply  of  a  Sigma/Timetal  21S 
(0/±45/90)s  laminate  subjected  to  HIP  and  isothermal  cyclic  loading 

Figure  34  Stress  average  computed  in  the  45°  ply  of  a  Sigma/Timetal  21S  (0/±45/90)s 
laminate  subjected  to  HIP  and  isothermal  cyclic  loading 

Figure  35  Fiber  stress  average  computed  in  the  45°  ply  of  a  Sigma/Timetal  21S 
(0/±45/90)s  laminate  subjected  to  HIP  and  isothermal  cyclic  loading 

Figure  36  Matrix  stress  average  computed  in  the  45°  ply  of  a  Sigma/Timetal  21S 
(0/±45/90)s  laminate  subjected  to  HIP  and  isothermal  cyclic  loading 

Figure  37  Stress  average  computed  in  the  90°  ply  of  a  Sigma/Timetal  21S  (0/±45/90)s 
laminate  subjected  to  HIP  and  isothermal  cyclic  loading 
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Figure  38  Fiber  stress  average  computed  in  the  90°  ply  of  a  Sigma/Timetal  21S 
(0/±45/90)s  laminate  subjected  to  HIP  and  isothermal  cyclic  loading. 

Figure  39  Matrix  stress  average  computed  in  the  90°  ply  of  a  Sigma/Timetal  21S 
(0/±45/90)s  laminate  subjected  to  HIP  and  isothermal  cyclic  loading 

Figure  40  Measured  and  computed  maximum  overall  axial  strains  in  Sigma/Timetal  21S 
laminates  under  isothermal  fatigue  loading  conditions 
shown  in  Table  2  for  the  slow  cycle  at  0.001  Hz 

Figure  41  Measured  and  computed  maximum  overall  axial  strains  in  Sigma/Timetal  21S 
laminates  under  isothermal  fatigue  loading  conditions 
shown  in  Table  2  for  the  fast  cycle  at  0.1  Hz 
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Figure  1.  Temperature  profile  obtained  for  a  Sigma/Timetal 

in  a  succeptor  furnace 
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Figure  2.  Applied  stress,  and  measured  total  strain  histories  in  a  multistep 

creep  test  for  a  Timetal  21 S  specimen  at  650°C 
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Figure  4.  Cyclic  loading  applied  to  a  Timetal  21S  specimen  at  650°C 
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Figure  5.  Maximum  and  minimum  strains  measured  under  the  cyclic  loading  of  Fig.  4 
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Maximum  and  minimum  strains  produced  under  each  of  the  cyclic  stress 

ranges  of  Fig.  4 
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Figure  7.  Tensile  stress-strain  data  measured  f< 

at  650°C  under  two  si 


Total  axial  strain  histories  measured  in  multistep  creep  tests  for  two 
(0)4  Sigma/Timetal  21S  laminates  at  650°C 
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Figure  10.  Total  axial  strain  history  measured  in  a  multistep  creep  test  for 
a  (0/±45/90)s  Sigma/Timetal  21S  laminate  at  650°C 
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Figure  12.  Total  axial  strain  history  measured  for  a  (0)4  Sigma/Timetal  21S  laminate 
at  650°C  under  the  fast  load  cycles  shown  in  inset 
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Figure  14.  Stress— axial  strain  response  measured  for  a  (0)4  Sigma/Timetal  21S 

laminate  under  the  fast  load  cycles  shown  in  inset  of  Fig.  12 
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Figure  18.  Stress— axial  strain  response  measured  for  a  (0/90)s  Sigma/Timetal  21S 

laminate  under  the  fast  load  cycles  shown  in  inset  of  Fig.  16 
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Figure  19.  Total  axial  strain  history  measured  for  a  (0/±45/90)s  Sigma/Timetal  21S 
laminate  at  650°C  under  the  slow  load  cycles  shown  in  inset 


Figure  20.  Total  axial  strain  history  measured  for  a  (0/±45/90)s  Sigma/Timetal  21S 
laminate  at  650°C  under  the  fast  load  cycles  shown  in  inset 
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Figure  23.  Maximum  stress— cycles  to  failure  (S— N)  data  measured  for 
Sigma/Timetal  21S  laminates  at  650°C  under  two  frequencies 


Figure  24.  Variation  of  secant  modulus  measured  in  two  (0)4  Sigma/Timetal  21S 

laminates  subjected  to  cyclic  loading  at  650°C 
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Figure  25.  Maximum  strains  measured  at  failure  in  fatigue  tests  of  Sigma/Timetal  21S 

laminates  at  650°C 
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Figure  28  Variation  of  the  transverse  normal  stress  in  fully  bonded  fibers  of 

a  Sigma/Timetal  21S  (0/90)s  laminate  with  the  overall  axial  stress 
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Figure  30  Isothermal  cyclic  loading  applied  to  a  Sigma/Timetal  21S  (0/±45/90)s  laminate 
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average  computed  in  the  45°  ply  of  a  Sigma/Timetal  21S  (0/±45/90) 
laminate  subjected  to  HIP  and  isothermal  cyclic  loading 
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Figure  35  Fiber  stress  average  computed  in  the  45°  ply  of  a  Sigma/Timetal  21S 
(0/±45/90)s  laminate  subjected  to  HIP  and  isothermal  cyclic  loading 
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Figure  38  Fiber  stress  average  computed  in  the  90°  ply  of  a  Sigma/Timetal  2 IS 
(0/±45/90)s  laminate  subjected  to  HIP  and  isothermal  cyclic  loading. 
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Figure  39  Matrix  stress  average  computed  in  the  90°  ply  of  a  Sigma/Timetal  21S 
(0/±45/90)s  laminate  subjected  to  HIP  and  isothermal  cyclic  loading 
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Figure  40  Measured  and  computed  maximum  overall  axial  strains  in  Sigma/Timetal  21S 
°  laminates  under  isothermal  fatigue  loading  conditions 

shown  in  Table  2  for  the  slow  cycle  at  0.001  Hz 
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